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1.0    SUMMARY 

IEC  International  Environmental  Consultants  Ltd.  was  retained  by  the  Ontario 
Ministry  of  the  Environment  to  operate  samplers  at  six  sites  in  the  Toronto 
area  during  the  28  day  province  wide  survey  of  PCB's  conducted  in  September 
-  October,  1979.  IEC  was  also  retained  to  conduct  data  analysis  for  all  25 
sites  sampled  in  the  survey  using  PCB  concentration  data  supplied  by  the 
Ministry  and  meteorologic  data  obtained  from  the  Ministry  and  Atmospheric 
Environment  Service  of  Environment  Canada. 

The  purpose  of  the  survey  was  to  provide  data  on  ambient  PCB  levels  at  the 
25  sites  which  were  selected  by  MOE  and  to  correlate  the  measured  PCB 
concentrations  with  meteorological  conditions  during  the  sampling  period. 

The  data  reduction  and  analysis  were  intended  to  demonstrate  whether  any 
significant  correlations  between  PCB  concentration  and  the  meteorological 
factors  wind  direction,  wind  speed,  and  air  temperature  existed. 

Linear  regression  correlations  between  PCB  concentration  and  meteorological 
parameters  were  performed.  Few  significant  correlations  between  PCB 
concentration  and  temperature  or  standard  deviation  of  temperature  were 
found.  Sixty-five  statistically  significant  correlations  between  PCB 
concentration  and  wind  direction  were  found,  while  thirty-four  significant 
correlations  between  PCB  concentration  and  wind  speed  were  obtained. 

Further  examination  of  eighteen  of  the  most  significant  correlations  showed 
that  while  all  were  statistically  significant,  the  data  were  skewed  such 

that  the  correlations  were  based  on  relatively  few  data  points. 
Consequently,  the  correlations  should  be  interpreted  as  trends. 
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2.0    INTRODUCTION 

In  June  1979  the  Ontario  Ministry  of  the  Environment  invited  proposals  for  a 
"Survey  of  Polychlon'nated  Biphenyls  in  Ambient  Air  in  the  Province  of 
Ontario".  The  intent  of  the  program  was  to  obtain  baseline  data  on  ambient 
PCB  concentrations  in  a  summer  survey  and  in  a  winter  survey  at  25  stations 
in  Ontario  and  to  evaluate  the  significance  of  correlations  between  PCB 
concentrations  and  meteorologic  data  at  each  station.  In  addition  to  this, 
there  was  to  be  an  intensive  survey  in  the  Mississauga  area  during  a  test 
burn  of  PCB-containing  oil  at  the  St.  Lawrence  Cement  Plant. 

The  program  was  subsequently  reduced  to  one  province-wide  survey  undertaken 
in  September-October,  1979.  This  report  presents  the  results  of  the 
evaluation  of  PCB  concentration  and  meteorologic  data  for  correlations, 
along  with  all  computer  output  for  the  September-October,  1979  survey.  Raw 
meteorologic  data  for  each  station  have  been  submitted  in  a  separate 
report. 
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3.0 


DATA  ACQUISITION  AND  REDUCTION 


3.1 


Meteorological  Data 


The  meteorological  data  used  in  the  regression  analyses  were  obtained  from 
the  Ontario  Ministry  of  the  Environment,  the  Atmospheric  Environment  Service 
of  Environment  Canada  and  Ontario  Hydro.  Meteorologic  data  were  obtained 
from  the  stations  listed  in  Table  1.  As  there  were  25  PCB  sampling  stations 
and  only  13  meteorologic  stations,  data  from  one  meteorologic  station 
frequently  was  used  at  more  than  one  PCB  sampling  station  for  purposes  of 
correlation  as  indicated  in  Table  2. 


TABLE  1:  METEOROLOGICAL  STATIONS 


Station  Name 

1. 

MOE-14016 

2. 

M0E-12032 

3. 

MOE-12034 

4, 

AES-LOND 

5. 

MOE-29026 

6. 

M0E-27049 

7. 

MOE-35033 

8. 

M0E-33003 

9. 

MOE-45025 

10. 

AES-KING 

11. 

M0E-77025 

12. 

AES-TBAY 

13. 

HYD-NANT 

Station  Address 

Opposite  Lambton  generating  station,  Hwy.  40,  Courtright 

Windsor  (not  used  -  temperature  data  not  available) 

C.P.  Tower,  Riverside  Drive,  Windsor 

London  Airport  (London-A) 

Brampton/Woodward  Avenue,  Hamilton 

Niagara  Falls  (Temp,  from  AES-St.  Catherines  "A") 

Evans/Arnold  Ave.  Etobicoke 

Lawrence/Kennedy  Ave.  Scarborough 

Ritson  Rd. /Olive  Ave.  Oshawa  (Temp,  from  AES,  Oshawa 

Water  Pollution  Control  Plant). 

Kingston  "A",  Kingston  Airport 

CKNC  TV  Tower,  Frood  Road,  Sudbury 

Thunder  Bay  "A",  Thunder  Bay  Airport 

Ontario  Hydro,  Nanticoke  Generating  Station 
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The  data  consisted  of  hourly  readings  of  temperature,  wind  speed  and  wind 
direction.  Where  meteorological  data  was  missing  at  a  particular  location 
the  missing  data  were  obtained  from  another  source.  Originally,  the  MOE 
station  at  Allanburg  was  to  be  used  rather  than  the  one  at  Niagara  Falls. 
However,  instrument  failures  resulted  in  an  unacceptably  high  percentage  of 
missing  data  (>50%).  Temperature  data  were  not  available  at  the  Niagara 
Falls  location,  but  were  obtained  from  the  AES  St.  Catherine's  "A"  Station. 
Similarly,  the  temperature  data  for  MOE-45025  (Oshawa)  were  obtained  from 
AES.  These  stations  provide  only  daily  temperature  maxima  and  minima, 
therefore  an  average,  calculated  for  each  day,   was  used. 

Wind  direction  data  are  coded  by  AES  into  36  classes  of  10  degrees  each. 
Thus  the  actual  degree  reading  is  approximately  equal  to  ten  times  the  class 
number.  Because  the  classes  used  by  AES,  and  by  IEC  in  the  p'reparation  of 
this  report,  were  different  (ie.  36  classes  of  10  degrees  by  AES  versus  16 
classes  of  22.5  degrees  by  IEC)  one  degree  was  subtracted  from  the  reading 
obtained  by  multiplying  AES  classification  by  10.  This  prevented  any 
intersection  of  class  limits  and  allowed  an  immediate  verification  that  wind 
direction  data  were  originally  coded. 

The  formula  used  to  translate  classes  back  to  degrees  was: 

Direction  =  (CODE  X  10)  -1 

This  gave  a  reading  in  degrees  which  was  approximately  centered  in  the  code 
class  defined  by  AES.  MOE  and  Ontario  Hydro  data  were  not  coded. 

Wind  speed  values  were  given  in  nautical  miles  per  hour  in  the  AES  Data.  A 
conversion  factor  of  1.852  (CRC,  1973)  was  used  to  convert  the  data  to  units 
of  km/hour.  MOE  and  Ontario  Hydro  data  were  available  in  metric  units. 

Copies  of  raw  meteorological  data  output  are  provided  in  a  separate  volume 
of  this  report. 
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TABLE  2:   LOCATION  OF  PCB  SAMPLING  SITES 


PCB  Station 


Location 


Met  Station 


BUR-S-U1 

Burlington,  Urban 

MOE  - 

29026 

HAM-S-I1 

Hamilton,  Industrial 

MOE  - 

29026 

HAM-S-I2 

Hamilton,  Industrial 

MOE  - 

29026 

HAM-S-U1 

Hamilton,  Urban 

MOE  - 

29026 

HAM-S-U2 

Hamilton,  Urban 

MOE  - 

29021 

KIN-S-U1 

Kingston,  Urban 

AES  - 

KING 

LON-S-U1 

London,  Urban 

AES  - 

LOND 

MIS-S-U1 

Mississauga,  Urban 

MOE  - 

35033 

MIS-S-U2 

Mississauga,  Urban 

MOE  - 

35033 

MIS-S-U3 

Mississauga,  Urban 

MOE  - 

35033 

MIS-S-U4 

Mississauga,  Urban 

MOE  - 

35033 

MOR-S-U1 

Mooretown,  Urban 

MOE  - 

12034 

NAN-S-R1 

Nanticoke,  Rural 

HYD  - 

NANT 

NAN-S-R2 

Nanticoke,  Rural 

HYD  - 

NANT 

0SH-S-U1 

Oshawa,  Urban 

MOE  - 

45025 

STC-S-R1 

St.  Catherines,  Rural 

MOE  - 

27049 

SAR-S-I1 

Sarnia,  Industrial 

MOE  - 

14016 

SAR-S-U1 

Sarnia,  Urban 

MOE  - 

14016 

SUD-S-S1 

Sudbury,  Suburban 

MOE  - 

77025 

SUD-S-R1 

Sudbury,  Rural 

MOE  - 

77025 

THU-S-R1 

Thunder  Bay,  Rural 

AES  - 

TBAY 

THU-S-U1 

Thunder  Bay,  Urban 

AES  - 

TBAY 

T0R-S-S1 

Toronto,  Suburban 

MOE  - 

33003 

TOR-S-U1 

Toronto,  Urban 

MOE  - 

33003 

WIN-S-U1 

Windsor,  Urban 

MOE  - 

12034 

Note:  U  =  Urban,  I  =  Industrial,  S  =  Suburban,  R  =  Rural 
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3.2  PCB  Data 

Blank-corrected  PCB  data  for  each  site  were  provided  to  IEC  by  MOE.  No 
modification  or  editing  was  performed.  PCB  station  names  and  locations  are 
given  in  Table  2. 

PCB  samplers  were  operated  on  a  24  hour  basis,  but  the  time  at  which 
sampling  started  and  stopped  varied  slightly  from  day  to  day  at  each  station 
and  varied  substantially  between  stations.  The  average  start  time  for  each 
station  was  determined  from  the  field  data  sheets  and  input  to  the  computer 
along  with  the  PCB  data.  This  was  done  to  achieve  as  close  a  match  between 
sampling  periods  (of  24  hours)  and  meteorological  data  as  possible. 

While  most  stations  started  the  study  period  on  the  27th  of  September,  two 
did  not.  These  were  BUR-S-U1,  which  started  on  the  26th  of  September,  and 
KIN-S-U1,  which  started  on  the  29th  of  September.  The  meteorological  data 
files  used  with  these  stations  were  modified  to  give  correct  results  with 
these  stations. 

Listings  of  the  PCB  data  files  used  are  provided  in  Appendix  1. 

3.3  Data  Reduction 

A  computer  program  was  developed  to  perform  the  data  analysis.   It  had 

provisions  for  entering  and  modifying  PCB  and  meteorological  data,  and 

performed  the  correlation  analysis  between  any  PCB  sampling  station  and  any 
meteorological  station.   A  complete  listing  of  the  program,  written  in 

BASIC,  is  provided  in  Appendix  2. 

Linear  correlation  analysis  was  used  to  calculate  correlation  coefficients. 
MOE  uses  similiar  programs  to  correlate  such  parameters  as  suspended 
particulate  matter  and  wind  direction. 
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3.4    Statistical  Calculations 

The  calculation  of  means,  standard  deviation,  linear  regression  constants 
and  correlation  coefficients  were  performed  using  standard  formulae.  They 
are   given  below: 

A.  Mean 

Mean  =  — ^ L 
n 

Xj  -  ith  data  point 

n  -  total  number  of  data  points 

B.  Standard  Deviation 
=  Axi  -  n(Mean)2 


Source  -  Barrowdale,  Roberts,  Ehle,  "Elementary  Computer 
Applications",  Toronto,  1974. 

C.  Correlation  Coefficient 

r  =    l  x^yi  -  n  EXj  Zy, 


(Ex-j2-  n  (Zxi)2)  X  (Eyi2-  n  (Eyi)2) 


Source  -  Barrowdale,  Roberts,  Ehle,  ibid. 
D.  Linear  Regression  Constants 


Slope  =  nExy  -  ExEy 
nEx2  -(Ex)2 


Source  -  Barrowdale,  Roberts,  Ehle,  ibid 
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The  t-test  was  used  to  evaluate  the  significance  of  the  correlation  coef- 
ficients calculated  for  each  data  set.  The  test  was  based  on  the  null 
hypothesis  that  the  calculated  regression  line  was  not  significantly  dif- 
ferent from  a  line  of  zero  slope  (representing  the  mean  value  of  all  "y" 
values) . 

The  formula  used  to  calculate  the  value  of  the  t-statistic  is  as  follows: 

t  =  rV  n  -  2 
1  -  r2 

Where:   r  =  correlation  coefficient 
n  =  number  of  data  points 

Source:  Volk,  "Applied  Statistics  for  Engineers",  Toronto,  1958. 

Thus,  if  a  correlation  is  found  to  have  t-value  greater  than  the  tabulated 
value  at  the  90*  confidence  limit  but  less  than  that  for  the  95%  confidence 
limit  then  the  null  hypothesis  is  false  at  the  90%  confidence  level,  but 
remains  true  at  the  95%  level  and  the  correlation  is  therefore  significant 
at  the  90%,  but  not  at  the  95%  level. 

There  are  two  possible  errors  which  may  result  from  any  statistical  test  of 
this  type,  known  as  Type  1  and  Type  2.  Type  1  are  those  in  which  the  null 
hypothesis  is  rejected  (ie  said  to  be  false)  when  it  is  not,  while  Type  2  is 
the  reverse,  that  is,  it  is  said  to  be  true  when  in  fact  it  is  false. 

Since  these  are  mirror  images  of  each  other,  it  follows  that  attempting  to 
eliminate  the  possibility  of  a  Type  1  error  may  result  in  a  Type  2  error. 
The  way  to  avoid  this  is  to  increase  the  sample  size. 

Statisticians  commonly  state  that  a  90%  significance  level  is  the  minimum 
level  for  a  correlation  to  be  accepted  as  significant,  while  the  95  and  99% 
levels  are  better,  and  tend  to  balance  Type  1  and  Type  2  errors  the  best,  as 
well.  Consequently,  these  three  levels  have  been  used  in  the  assessment  of 
the  significance  of  the  correlations  between  PCB  concentrations  and  meteor- 
ologic  data. 


The  correlations  to  be  performed  were: 

1)  PCB  concentration  vs  wind  direction  class 

2)  PCB  concentration  vs  wind  speed 

3)  PCB  concentration  vs  ambient  temperature 

4)  PCB  concentration  vs  standard  deviation  of  daily  temperature 

The  approach  used  in  the  program  was  as  follows: 

1)  load  PCB  data  file  for  one  station  and  sampling  start  time 

2)  load  corresponding  meteorological  data  file 

3)  calculate  daily  average  temperatures,  with  "day" 
starting  at  start  time  for  the  PCB  station 

4)  calculate  standard  deviation  of  daily  temperatures  from  daily  mean 
temperatures 

5)  classify  wind  direction  into  16  classes 

6)  classify  wind  speeds  into  five  classes 

7)  normalize  absolute  counts  (number  of  hours  per  wind  direction) 
obtained  in  (5)  and  (6)  above  to  %   time 

8)  perform  linear  regression  analysis  as  outlined  above 

9)  output  results  in  tabular  form 

10)  output  wind  direction  results  in  histogram  form. 

Output  from  all  stations  is  provided  in  Appendix  3  of  this  report.  The 
program  requires  approximately  40  Kbytes  of  memory  to  run,  divided  between 
program  storage,  17  Kbytes  and  data  storage,  23  Kbytes. 
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I  EC 

4.0  RESULTS  AND  DISCUSSION 

4.1  PCB  Data 

The  means,  standard  deviations  and  maximum  and  minimum  readings  for  the  raw 
PCB  data  provided  by  MOE  for  all  stations  are  shown  in  Table  3.  The  highest 
one  day  PCB  concentration  of  74  ng/m-3  was  found  in  Hamilton.  The 
highest  average  reading  of  9.8  ng/m^  was  obtained  in  downtown  Toronto. 
The  lowest  one-day  concentration  of  36  ng/m^  was  found  in  Kingston,  while 
the  lowest  average  of  3.8  ng/m^  was  found  at  the  rural  site  in  the 
Thunder  Bay  area.  All  the  measured  daily  concentrations  are  below  the 
Ministry  guidelines  of  150  ng/rn^. 

In  general,  the  standard  deviation  for  each  station  was  high,  indicating  a 
large  day-to-day  variability  in  ambient  PCB  concentration.  While  there  was 
only  one  station  in  which  standard  deviation  exceeded  the  mean,  in  most 
cases  the  standard  deviation  was  at  least  66%  of  the  mean. 

A  map  of  Ontario  indicating  the  25  sampling  sites  is  presented  in  Figure  1. 
The  mean  and  standard  deviation  of  the  PCB  concentrations  at  each  site 
during  the  study  are  shown.  Daily  concentrations  for  all  stations  are  shown 
on  similiar  maps  in  Appendix  4. 

4.2  Correlations 

The  use  of  a  simple  linear  regression  analysis  to  determine  possible 
relationships  between  the  PCB  concentration  and  the  various  meteorological 

parameters  is  not  a  statistically  rigorous  treatment  of  the  data.  To  be 
accurate,  multiple  correlation  analysis  incorporating  all  of  the 
meteorological  parameters  in  the  study  would  have  to  be  performed.  This 
would  ensure  that  there  are  no  internal  correlations,  for  instance,  between 
wind  direction  and  wind  speed.  If  this  were  the  case,  then  a  strictly 
rigorous  treatment  of  the  problem  would  demand  that  the  statistically  less 
significant  variable  be  discarded  and  the  regression  re-run.  Meteorologists 
commonly  assume  that  there  is  no  significant  relation  between*  the  various 
parameters  on  a  large  scale  over  a  long  period  of  time.   A  28  day  time 
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TABLE  3: 

RAW  PCB  DATA  STATISTICS 

STATION 

#  VALID 
OBSERVATIONS 

MIN. 
(ng/mJ) 

MAX  ~ 
(ng/mJ) 

MEAN. 
(ng/mJ) 

STANDARD 
DEVIATION 

BUR-S-UI 

25 

1.0 

24 

5.6 

5.2 

HAM-S-11 

26 

1.1 

22 

6.9 

5.6 

HAM- S- 12 

25 

1.0 

18 

5.9 

4.6 

HAM-S-UI 

25 

1.0 

22 

5.7 

5.0 

HAM-S-U2 

26 

1.2 

74 

9.2 

13.9 

KIN-S-UI 

24 

0.4 

10 

4.7 

3.0 

LON-S-UI 

18 

0.4 

18 

8.1 

4.8 

MIS-S-UI 

24 

0.9 

15 

5.3 

3.9 

MIS-S-U2 

25 

0.8 

15 

4.9 

3.9 

MIS-S-U3 

25 

0.6 

21 

6.3 

5.7 

MIS-S-U4 

27 

0.9 

15 

5.3 

3.6 

MOR-S-U1 

26 

0.6 

16 

4.8 

3.6 

NAN-S-R1 

27 

1.0 

14 

5.2 

3.7 

NAN-S-R2 

27 

0.9 

18 

5.1 

4.3 

OSH-S-U1 

26 

1.0 

14 

5.7 

4.4 

STC-S-RI 

26 

1.6 

23 

6.6 

5.6 

SAR-S-II 

24 

0.6 

21 

7.3 

5.6 

SAR-S-UI 

24 

1.8 

47 

9.8 

9.4 

SUD-S-SI 

24 

1.0 

12 

4.4 

2.9 

SUD-S-R1 

26 

1.7 

17 

6.8 

4.7 

THU-S-R1 

26 

0.7 

10 

3.8 

2.6 

THU-S-UI 

27 

1.3 

13 

5.7 

4.2 

TOR-S-SI 

23 

1.6 

15 

5.3 

3.3 

TOR-S-UI 

25 

1.9 

32 

9.6 

7.3 

WIN-S-UI 

27 

1.5 

16 

5.6 

3.6 
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period  should  minimize  these  problems  as  well  as  provide  an  adequate  data 
base  to  begin  to  define  ambient  and  background  levels  of  PCB's  in  the 
Province  of  Ontario. 

4.2.1  Correlations  Between  PCB  Concentration  and  Temperature 

A  total  of  five  significant  correlations  with  average  daily  temperature  were 

found;   only  two  of  these  were  significant  at  the  95%  level,  while  the 

others  were  significant  at  the  90%  level.  Table  4  summarized  all  the 
calculated  correlation  coefficients. 

Correlations  with  standard  deviation  of  daily  temperature  were  performed  to 
provide  a  rough  estimate  of  the  behaviour  of  PCB's  contained  in  the  atmo- 
sphere under  varying  conditions  of  atmospheric  stability.  It  was  antici- 
pated that  few  significant  correlations  would  be  found  and  that  those  that 
did  occur  would  have  negative  correlation  coefficients.  Increasing  standard 
deviation  of  temperature  would  indicate  less  atmospheric  stability,  and  thus 
greater  mixing  between  boundary  layer  and  upper  atmosphere,  with  resultant 
dilution  of  PCB  concentration.  The  calculated  correlation  coefficient  are 
given  in  Table  4.  Of  the  three  significant  correlations  found  only  one  had 
the  expected  negative  correlation  coefficient. 

4.2.2  Correlations  Between  PCB  Concentration  and  Wind  Direction  and  Speed 

The  calculated  correlation  coefficients  for  PCB  concentration  with  wind 

direction  classes  are  presented  in  Table  5a.  Those  correlations  which  are 
statistically  significant  are  presented  in  Table  5b.   A  summary  of  the 

number  of  statistically  significant  correlations  by  wind  direction  and  level 
of  significance  is  presented  in  Table  6. 

The  values  of  the  t  statistic  necessary  for  the  correlations  to  be 
significant  at  the  90%,  95%  and  99%  levels  are  given  in  Table  7.  In  this 
case,  a  significance  of  90%  means  that  in  one  out  of  10  cases,  a  correlation 
rejected  as  insignificant  will  actually  be  significant.  (ie.  a  Type  2 
error). 
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TABLE  4: 


CORRELATION  COEFFICIENTS  -  PCB  CONCENTRATION  vs  TEMPERATURE 


STATION 


vs  AVERAGE 
TEMPERATURE 


vs  STD.   DEV. 
TEMPERATURE 


BUR-S-U1 
HAM-S-I1 
HAM-S-I2 
HAM-S-U1 
HAM-S-U2 
KIN-S-U1 
L0N-S-U1 
MIS-S-U1 
MIS-S-U2 
MIS-S-U3 
MIS-S-U4 
MOR-S-U1 
NAN-S-R1 
NAN-S-R2 
OSH-S-U1 
STC-S-R1 
SAR-S-I1 
SAR-S-U1 
SUD-S-S1 
SUD-S-R1 
THU-S-R1 
THU-S-U1 
TOR-S-S1 
TOR-S-U1 
WIN-S-U1 


0.0515 

0.0662 

0.2656 

0.2065 

0.3026 

0.0604 

0.4046 

0.3772" 

0.3819 


0.4115 

0.0561 

0.0620 

0.2734 

0.2593 

0.3630 

0.1142 

0.45742 

0.1983 

0.0247 

0.2400 

0.0257 

0.  36431 

0.3391 

0.3218 

0.1847 


-0.0711 

0.0593 

2.1025 

0.0185 

-0.1843 

-0.0784 

-0.1239 

0.0824 

0.1799 

-0.0368 

0.4323' 

-0.1325 

0.0477 

0.3069 

0.0234 

-0.4103' 

-0.0901 

0.0722^ 

0.4341' 

-0.0089 

-0.0043 

0.1211 

0.0248 

-0.1867 

0.0021 


NOTES  -   1     Significant   at  90%  level 
2     Significant   at  95%  level 
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TABLE  5a:   CORRELATION  COEFFICIENTS  -  PCB  CONCENTRATION  vs  WIND  DIRECTION  CLASS 


STATION 


UNU 


NW 


WNW 


5  5. 


BUR-S-U1 

0.2958 

-0.37121 

-0.1369 

C.4324 

3 
0.6227 

3 

0.8443 

-0.0171 
0.3931 

-0.0713 

HAK-S-I1 

0.3160 

-0.0727 

-0.0482 

0.2605 
0.5093 

0.1S95 

3 
0.5299 

0.0785 

0.0727 

HAM-S-I2 

C.7321  3 

-0.0010 

0.1681 

0.1426 

-0.0fc20 

-0.1615 

HAM-S-U1 

0.72273 

-0.1344 

0.1613 

0.45992 

0.5436 

0.3327 

0.1250 

-O.IX.43 

HAT-S-U2 

-0.0935 

-0.0636 

-0.1543 

-0.0139 

-0.0055 

-0.0951 

-0.0977 

-0.0ol7 

KI N-S-U 1 

-0.2596 

-C.I219 

0.1512 

-0.2005 

0.1367 

-0.0708 

0.0359 

0. 12ol 

lON-S-Ul 

0.3625 

0.2424 

0.2488 

C.3447 

0.1954 

0.3767 

-0.0370 

-0.0100 

MIS-S-U1 

0.2094 

-0.0125 

0.0424 

0.0008 

-0.294S 

0.01^2 

0.1261 

0.73Jo3 

MIS-S-U2 

0.1722 

-0.1291 

-0.1151 

-0.0338 

-0.3041 

-0.1594 

-0.1238 

0.52963 

K1S-S-U3 

0.2654 

0.1986 

0.2539 

0.0018 

-0.3110 

-0.1617 

0. 108o 

0.  5o92 

M1S-S-U4 

0.0261 

-0.2145 

-0.2132 

-0.0110 

-0.0629 

-0.046o 

-0.0ot>2 

0.2137 

MOR-S-Ul 

0.0074 

-0.1259 

C.0522 

-0.0446 

0.2403 

-0.0706 

0.0871 

0.63123 

NAN-S-R1 

0.2930 

0.1346 

0.1011 

0.0714 

-0.0267 

-0.0404 

-0.0903 

-0.1094 

NAN-S-R2 

0.1686 

0.2673 

0.42452 

0.36601 

0.1329 

0.0539 

0.0779 

0.0091 

OSH-S-U1 

-0.0339 

0.1694 

0.1031 

0.2262 

0.3041 

-0.1369 

-0.0831 

-0.0419 

STC-S-R1 

0.46712 

0.2050 

0.0226 

0.0108 

-0.2817 

0.39282 

0.42o82 

0.1035 

SAR-S-I1 

0.1560 

0.63923 

0.55413 

0.1033 

0.1906 

-0.1124 

-0.0809 

0. 14o0 

SAR-S-U1 

0.1910 

0.68943 

0.48272 

-C.0768 

0.0127 

•0.1419 

-0.2020 

-0.C449 

SUD-S-S1 

0.3011 

0.2631 

0.36591 

0.2370 

0.0455 

0.2379 

0.0063 

-0.2234 

SUO-S-R1 

-0.1342 

-0.0731 

-0.0098 

-0.0150 

-0.2319 

-0.0415 

-0.2562 

-0.2153 

THU-S-R1 

0.0403 

-0.0242 

-0.3871 ] 

-0.2591 

-0.5231 3 

-0.41272 

-0.2781 

0.0o76 

THU-S-U1 

0.3096 

-0.0784 

-0.44692 

-0.34  961 

-0.46452 

-0.2287 

-0.0330 

0.?Ot»Z 

TOR-S-S1 

0.5794  3 

C.42802 

0.36821 

-0.0724 

-0.0075 

0.52052 

C.  1 704 

0.2928 

TOR-S-U1 

0.1983 

0.1952 

0.48562 

0.43772 

0.1238 

0.35091 

-0.0417 

0.43O02 

W1N-S-U1 

0.0423 

0.0026 

0.38282 

0.3436  ] 

0.63043 

0.1538 

0.3062 

0.2074 

STATION 


SSE 


SE 


USE 


ENE 


m 


NM 


BUR-S-U1 

0.0287 

-0.2513 

-0.1657 

-0.0449 

0.2071 

0.0612 

-0.0181 

-0.0791 

HAM-S-11 

0.0844 

-0.1929 

-0.0904 

0.1598 

-0.1896 

-0.1388 

-0.0348 

-0.1115 

HAM-5-I2 

0.0728 

-0.1704 

-0.3172 

0.0062 

-0.3012 

-0.1304 

0.19O6 

0.0726 

HAM-S-U1 

0.1082 

-0.2781 

-0.36501 

-0.0412 

-0.2122 

-0.1047 

0.1227 

0.04  70 

HAM-S-U2 

-0.1304 

0.53643 

0.2755 

-0.1530 

-0.1837 

-0.1203 

-0.0167 

-O.OOoS 

KIN-S-U1 

0.1135 

0.0436 

0.0778 

-0.1389 

-0.2572 

-0.1907 

0.2087 

-0.21o2 

L0N-S-U1 

0.1689 

0.0438 

-0.2064 

-0.0770 

-0.2976 

-0.3088 

0.0429 

-0.0640 

M1S-S-U1 

0.51332 

-0.04  55 

-0.2863 

-0.1776 

-0.0137 

-0.0471 

-0.0998 

0.0728 

MIS-S-U2 

0.6121 3 

-0.0003 

-0.1503 

-0.1216 

0.0380 

0.2083 

-0.0496 

0.0023 

MIS-S-U3 

0.55173 

-0.0146 

-0.2070 

-0.2442 

0.0686 

0.0133 

-0.1703 

0.0765 

MIS-S-U4 

0.38432 
0.45252 

-0.0349 

-0.1555 

-0.2299 

0.0039 

0.0881 

0.1416 

0.1702 

M0R-S-U1 

-0.2678 

-0.1002 

0.0049 

-0.2279 

-0.1216 

-0.0854 

-O.U006 

NAN-S-R1 

0.0519 

-0.1448 

-0.1309 

0.0047 

-0.0026 

-0.0529 

-0.0689 

0.1810 

NAN-S-R2 

0.1301 

-0.0829 

-0.3787  * 

-0.1546 

-0.2362 

-0.2168 

0.0602 

0.19U 

0SH-S-U1 

-0.3666  1 

0.0138 

0.4367  2 

0.0741 

0.1455 

-0.0038 

-0.0182 

0.06o5 

5TC-S-R1 

0.0291 

0.0174 

-0.2176 

0.0128 

-0.2502 

-0.0728 

0.0045 

-0.2394 

SAR-S-I1 

0.0745 

-0.2743 

-0.2042 

-0.1562 

-0.2680 

0.0065 

-0.1695 

0.1347 

SAR-S-U1 

0.1293 

-0.3493  J 

-0.1657 

-0.1234 

-0.1982 

-0.1359 

-0.0753 

-0.0309 

0.3026 

SUD-S-S1 

-0.2118 
0.43122 
0.3700  1 

-0.1939 

-0.0419 

-0.1029 

-0.1478 

0.0441 

0.07U6 

SU0-S-R1 

0.1761 
0.4745 2 
0.4027 

0.2421 

0.0929 

0.2845 

0.36441 

-0.0758 

O.OOod 

THU-S-R1 

C.41682 

0.3393  l 

0.1910 

0.52843 

0.1518 

0.3323 

THU-S-Ul 

0.2891 

0.2172 

0.2400 

0.0811 

0.3I641 

0.3&492 

0.2421 

T0R-S-S1 

0.1021 

-0.1051 

-0.3156 

-0.43612 

-0.1968 

-C.1M8 

-0.1761 

0.O504 

T0R-S-U1 

-0.0578 

-0.2699 

-0.38051 

-0.2770 

-0.0070 

-0.1297 

-0.2752 

O.U204 

WIN-S-U1 

0.0064 

-O.2e50 

-0.2889 

-o.ieio 

-0.0758 

-0.0729 

0.0360 

-0.0>43 

Notes:  1  Significant  at  901  level 
'  Significant  at  95J  level 
3     Significant   at   99;   level 
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TABLE  5b:   SIGNIFICANT  CORRELATIONS  -  PCB  CONCENTRATIONS  vs  WIND  DIRECTION 


STATION 


Im 


WNW 


WSW 


sw 


ss- 


BUR-S-U1 
HAM-S-I1 
HAM-S-I2 
HAM-S-U1 
HAM-S-U2 
K1N-S-U1 
L0N-S-U1 
M1S-S-U1 
M1S-S-U2 
MIS-S-U3 
MIS-S-U4 
MOR-S-Ul 
NAN-S-R1 
NAN-S-R2 
OSH-S-U1 
STC-S-R1 
SAR-S-U 
SAR-S-U1 
SUD-S-S1 
SUD-S-R1 
THU-S-R1 
THU-S-U1 
TOR-S-S1 
TOR-S-U1 
WIN-S-U1 


0.7321 
0.7227 


0.4671 


0.5794 


-0.3712 


0.6392 
0.6694 


0.425 


0.4245' 


0.5541 J 
0.48272 
0.36591 


-0.3871 ] 

-0.44692 

0.36821 

Q.48S62 

0.38282 


1.4324' 


0.5293' 
0.4599' 


0.3560 


-0.3496 


0-43772 

0.34361 


73 


0.6227 


0.5299 
0.5436 


-0.5231 

-0.4645' 


0.6304' 


0.8448 


0.3928' 


-0.41272 
-2 


0.5205' 
0.3509 


0.3931' 


0.42682 


0.73063 
Q.62933 
0.56923 

0.63123 


0.4360 


STATION 


SSE 


SE 


ESE 


EN£ 


NE 


NNE 


BUR-S-U1 
HAM-S-I1 
MAM-S-12 
HAM-S-U1 
HAM-S-U2 
KIN-S-U1 
L0N-S-U1 
MIS-S-U1 
MIS-S-U2 
MIS-S-U3 
K1S-S-U4 
MOR-S-Ul 
NAN-S-R1 
NAN-S-R2 
0SH-S-U1 
STC-S-R1 
SAR-S-I1 
SAR-S-U1 
SU0-S-S1 
SU0-S-R1 
THU-S-R1 
THU-S-Ul 
TOR -S -SI 
T0R-S-U1 
WIN-S-U1 


0.51332 
0.61213 
0.55173 
0.38432 
0.45252 


-0.3666 


-0.3493 


-0.4745 
0.4027 


0.53643 


0.4312 
0.3700 


-0.3650 


-0.3787 
0.4367 


0.4168 


-0.3805 


0.3393 


-0.4361 


0.  3844  > 
0.5284 


0.3364 


0.33231 


0.3849' 


Notes:      1     Significant  at  90S  level 

?     Significant   at  951  level 
3     Significant   at  991  level 
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TABLE  6:   SIGNIFICANT  CORRELATIONS  -  PCB  CONCENTRATION  vs  WIND  DIRECTION 


SIGNIFICANCE  WIND  DIRECTION 

(%)     N    NNW    NW   WNW    W    WSW    SW    SSW    S    SSE    SE    ESE    E    ENE   NE    NNE   TOTAL 

90     -     1     33-1     -_2     131-     2-1     18 

95      1      1      5     3    13      2-5      1     2     1     -  1     -     29 

18 

1     1     65 


99 

3 

2 

1 

1 

5 

1 

- 

4 

2 

1 

- 

- 

- 

1 

- 
TOTAL 

4 

4 

9 

7 

6 

5 

2 

4 

9 

3 

5 

2 

0 

3 

TABLE  7:    ABSOLUTE  VALUE  OF  t-STATISTIC  NECESSARY  FOR 
GIVEN  SIGNIFICANCE  LEVELS 


#  VALID 

90% 

95% 

99% 

DAYS 

SIGNIFICANCE 

SIGNIFICANCE 

SIGNIFICANCE 

20 

1.734 

2.181 

2.878 

21 

1.729 

2.593 

2.861 

22 

1.725 

2.086 

2.845 

23 

1.721 

2.888 

2.631 

24 

1.717 

2.074 

2.819 

25 

1.714 

2.369 

2.807 

26 

1.711 

2.064 

2.797 

27 

1.708 

2.060 

2.787 

28 

1.706 

2.056 

2.779 
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Table  8  presents  the  correlation  coefficients  found  for  PCB  concentration 
versus  wind  speed  and  Table  9  presents  a  survey  of  the  statistically 
significant  correlations. 

When  interpreting  the  results  of  correlation  analyses,  care  must  be  taken  to 
avoid  attributing  greater  importance  to  the  correlations  than  is  justified. 
As  an  example,  at  station  SAR-S-U1  the  correlation  coefficient  for  PCB 
concentration  versus  wind  from  the  north-northwest  is  0.6894,  and  is 
significant  at  the  99%  level.  Examination  of  a  graph  of  the  data,  included 
as  Figure  2,  shows  that  the  majority  of  the  points,  are  on  the  ordinate  of 
the  graph.  Thus  the  data  are  skewed  towards  the  lower  end  of  the  graph,  and 
the  correlation  is  being  based  primarily  on  the  four  other  points  on  the 
graph.  Consequently,  the  correlations  should  be  interpreted  as  trends  in 
most  cases.  The  solution  to  this  type  of  problem  is  a  larger  data  base. 
Plots  of  a  representative  number  of  the  most  significant  correlations  are 
given  in  Appendix  5. 

In  general,  sources  of  PCB's  may  be  characterized  as  local,  medium  range  and 
long  range.  To  identify  a  local  source  conclusively,  it  must  be  encompassed 
by  sampling  stations,  thus  allowing  observation  during  all  wind  conditions. 
The  correlation  with  wind  speed,  however,  can  provide  an  indication  of  the 
distance  of  the  sample  site  from  the  source.  A  significant  positive 
correlation  between  PCB  concentration  and  a  single  wind  direction,  combined 
with  a  significant  positive  correlation  with  wind  speeds  of  0-5  km/hr  may 
imply  local  or  medium-range  sources  of  PCB's.  As  internal  correlations  may 
exist  between  wind  direction  and  speed,  it  is  not  possible  to  draw  more 
definite  conclusions  without  undertaking  multiple  regression  or  vector 
analysis  of  the  data.  This  was  beyond  the  scope  of  this  project. 

No  province-wide  trends  in  correlations  of  PCB  concentrations  with  transport 
of  PCB's  based  on  the  data  presented  in  this  report  because  of  potential 
interferences  from  local  meteorologic  conditions.  As  discussed  in  this 
section,  however,  there  are  definite  trends  suggesting  short  and  medium  range 
transport  at  many  stations.  The  use  of  the  histogram  plots  produced  by  the 
computer  program  makes  identification  of  these  trends  relatively  easy.  A 
discussion  of  trends  at  each  sampling  station  follows. 
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TABLE   8:      CORRELATION  COEFFICIENT   -   PCB  CONCENTRATION   vs  WIND  SPEED 


WIND  SPEED   (km/hr.) 


STATION 


0-5 


5-10 


10-15 


15-20 


>20 


BUR-S-U1 

0.2825 

0.0596 

0.0241 

-0.874 

-0.2678 

HAM-S-I1 

0.3116 

0.1667 

0.2997 

-0.2075 

-0.46172 

HAM-S-I2 

0.62043 

0.0674 

0.6818 

-0.38981 

-0.39451 

HAM-S-U1 

0.68983 

-0.0453 

-0.0038 

-0.42312 

-0.37921 

HAM-S-U2 

-0.0453 

-0.2067 

-0.42222 

-0.2502 

0.5517 

KIN-S-U1 

0.39781 

-0.0679 

-0.39951 

-0.2705 

0.0970 

L0N-S-U1 

0.0740 

0.3181 

0.1910 

-0.0953 

-0.3011 

MIS-S-U1 

0. 2650 

-0.1669 

-0.0348 

-0.2510 

0.0317 

MIS-S-U2 

0.0648 

-0.0576 

0.0586 

-0.1168 

0.0129 

MIS-S-U3 

0.2094 

0.0517 

-0.1781 

-0.3142 

0.0397 

MIS-S-U4 

-0.1753 

-0.0267 

0.1475 

-0.0303 

0.51893 

M0R-S-U1 

0.1761 

-0.0703 

0.020 7 

-0.1090 

-0.0907 

1 

NAN-S-R1 

0.4108 2 

0.1790 

-0.055b 

-0.2350 

-0.3278 

NAN-S-R2 

0.3120 

0.33831 

-0.1065 

-0.2508 

-0.34081 

0SH-S-U1 

0.1502 

-0.1700 

-0.0690 

-0.0289 

-0.0498 

STC-S-R1 

0.45582 

-0.38942 

-0.33831 

-0.2015 

-0.0820 

SAR-S-I1 

0.1303 

0.2478 

0.2160 

-0.1676 

-0.2769 

SAR-S-U1 

-0.0028 

0.1189 

0.43672 

-0.0639 

-0.2548 

SUD-S-S1 

-0.2911 

0.2502 

0.2360 

0.1466 

-0.1300 

SUD-S-R1 

-0.36251 

-0.40822 

-0.35321 

0.2445 

0.53203 

THU-S-R1 

0.43282 

0.2115 

-0.39002 

-0.35831 

-0.2125 

THU-S-U1 

0.47882 

-0.0227 

-0.41412 

-0.2886 

-0.0181 

T0R-S-S1 

0.3130 

-0.0217 

-0.1559 

-0.2261 

-0.1976 

T0R-S-U1 

0.51773 

-0.0177 

-0.39751 

-0.3325 

-0.3139 

WIN-S-U1 

0.57753 

0.1731 

-0.47692 

-0.41942 

-0.2981 

Notes:      *     Significant   at  90%  level 

2  Significant  at  95%  level 

3  Significant   at  99%  level 
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TABLE    9     SIGNIFICANT  CORRELATIONS   -  PCB  CONCENTRATION  vs  WIND  SPEED 


SIGNIFICANCE 

(%)     0-5 


WIND  SPEED  (km/hr.) 
5-10     10-15     15-20     >20      TOTAL 


90 


95 


4       13 


99 


TOTAL 


3        9 


4        8 


7 
34 
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BUR-S-U1 

A  definite  trend  is  observed  when  winds  are  coming  from  the  west  towards  the 
station.  The  largest  correlation  coefficient  observed,  0.8448,  was  found  at 
this  station.  There  is  a  large  electrical  substation  to  the  south-southwest 
of  the  station.  As  the  Niagara  escarpment  is  northwest  of  the  station,  it 
is  possible  that  this  resulted  in  local  wind  direction  being  different  from 
that  recorded  at  Hamilton  airport  from  which  the  meteorological  data  was 
obtained.  Graphs  of  the  data  for  the  west  and  west-southwest  directions  are 
provided  in  Appendix  5.  There  were  no  significant  correlations  with  wind 
speed. 

HAM-S-I1 

A  95%  signficance  level  correlation  was  found  with  winds  from  the  southwest. 
A  general  trend  may  also  be  inferred  with  winds  ranging  from  WNW  through  S. 
A  95%  level  correlation  was  found  with  winds  in  the  >20  km/hr  class.  As 
this  was  a  negative  correlation,  it  implies  that  cleaner  air  was  coming  into 
the  area  during  periods  of  high  winds.  This  implies  accumulation  of  PCB's 
from  anthropogenic  sources  in  the  urban  envelope,  with  dispersal  during 
periods  of  high  winds. 

HAM-S-I2 

A  high  correlation  coefficient  was  found  when  winds  were  from  the  north.  A 
trend  was  also  found  from  the  NW  to  the  WSW  directions.  The  largest 
correlation  coefficient  in  that  quadrant  was  with  winds  from  the  west.  When 

winds  were  from  the  east,  a  negative  correlation  coefficient  was  found. 
This  suggests  local  or  medium  range  sources.  This  is  supported  by  a  99% 
level  correlation  with  wind  speeds  in  the  0-5  km/hr  class.  It  is  further 
supported  by  the  90%  level  negative  correlations  with  wind  speeds  in  the 
15-20  and  >20  km/hr  classes,  which  suggests  that  cleaner  air  enters  the  area 
during  periods  of  high  winds.  A  graph  of  the  data  for  wind  from  the  north 
and  the  0-5  km/hr  wind  speed  class  is  provided  in  Appendix  5. 
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HAM-S-U1 

The  same  trends  and  patterns  as  were  found  at  HAM-S-I1  were  found  at  this 
station.  A  plot  of  data  for  wind  from  the  north  and  the  0-5  km/hr  wind 
speed  class  is  given  in  Appendix  5. 

HAM-S-U2 

In  contrast  to  the  previous  stations,  only  the  south-southeast  direction 
showed  a  significant  correlation,  which  may  indicate  local  sources.  The  >20 
km/hr  class  showed  a  99%  level  positive  correlation,  suggesting  that  this 
station  may  be  further  from  the  source  than  the  previous  stations,  as  well. 
A  plot  of  the  >20  km/hr  wind  speed  class  data  is  given  in  Appendix  5. 

KIN-S-U1 


This  is  an  excellent  example  of  a  substantially  clean  area.  There  is  no 
significant  pattern  apparent  from  examination  of  the  histogram,  but  the  0-5 
and  10-15  km/hr  classes  both  showed  correlations  at  the  90%  level.  This 
could  be  due  to  the  skewed  data  problem  outlined  earlier  or  it  could  be 
local  sources.  The  west-east  pair  show  that  a  source  may  be  located  to  the 
west  of  the  site,  but  the  correlations  are  not  statistically  significant, 
and  so  no  definite  conclusion  can  be  drawn  from  the  data. 

L0N-S-U1 

While  no  significant  correlations  were  found  at  this  station,  the  pattern 
suggests  that  PCB's  may  be  transported  from  the  Sarnia  area.  This  is 
supported  by  a  corresponding  pattern  at  the  Sarnia  sites. 

MIS-S-U1,  MIS-S-U2,  MIS-S-U3,  MIS-S-U4 

A  trend  with  winds  ranging  from  S  to  SW  is  evident  at  these  stations,  with 
the  south-southwest  direction  showing  the  highest  correlation  coefficient  of 
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0.7306.  All  of  the  stations  in  Mississauga  showed  this  pattern.  This  may 
be  indicative  of  a  local  source  generally  to  the  south-southwest  of  the 
stations,  as  the  stations  in  Toronto  (T0R-S-U1  and  T0R-S-S1)  did  not  show 
this  pattern  under  similar  meteorological  conditions.  A  report  prepared  by 
the  Ontario  Research  Foundation  for  the  Ministry  entitled  "An  Assessment  of 
Atmospheric  PCB  Levels  in  Mississauga"  found  similar  results  using  sampling 
sites  yery  near  to  the  ones  used  in  this  study.  Plots  of  the  data  for  winds 
from  the  south  at  stations  -U2  and  -U3,  and  from  the  south-southwest  at 
stations  -Ul  and  -U3  are  given  in  Appendix  5.  Only  MIS-S-U4  had  a 
significant  correlation  with  wind  speed,  this  being  in  the  >20  km/hr  class 
at  a  significance  level  of  99%. 

M0R-S-U1 

A  trend  indicating  sources  to  the  south-southwest  was  found  at  this  station, 
with  the  south-southwest  correlation  having  a  99%  significance  level.  No 
significant  correlations  with  wind  speed  were  found. 

NAN-S-R1 

While  there  were  no  significant  correlations  with  wind  direction  at  this 
station,  a  trend  was  observed  with  winds  coming  from  the  NNE  through  to  WNW. 
A  95%  significance  level  correlation  was  found  for  the  0-5  km/hr  wind  speed 
class  and  a  90%  level  negative  correlation  was  found  for  the  >20  km/hr 
class.  These  trends  suggest  the  possibility  of  medium  range  transport  from 
a  generally  northerly  direction. 

NAN-S-R2 


At  this  station,  the  pattern  observed  at  NAN-S-R1  is  repeated,  but  the 
correlations  with  winds  from  the  northeast  and  west-northwest  had  95%  and 
90%  significance  levels,  respectively.  The  south-southeast  showed  a  90% 
level  negative  correlation.  This  tends  to  support  the  trend  inferred  above 
for  the  northerly  direction  and  indicates  relatively  cleaner  air  enters  the 
area  from  the  southeast. 
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• 


0SH-S-U1 

There  was  a  95%  level  correlation  with  wind  direction  from  the  southeast  at 
this  station,  suggesting  local  sources.  A  trend  was  also  apparent  in  the 
NNW  to  W  directions,  although  the  correlations  were  not  statistically 
significant. 

STC-S-R1 

The  correlation  data  indicate  local  sources  to  the  north  and  possibly  to  the 
west-southwest  and  south-southwest  of  this  station.  All  three  directions 
had  correlations  significant  at  the  95%  level.  The  0-5  km/hr  wind  speed 
class  showed  a  95%  significance  level  correlation  while  the  5-10  and  10-20 
km/hr  classes  showed  95%  and  90%  level  significant  correlations, 
respectively,  indicating  short  to  medium  range  transport. 

SAR-S-I1,  SAR-S-U1 

Both  of  these  stations  showed  trends  indicating  PCB  sources  to  the 
north-northwest.  The  generally  east-southeast  direction  showed  negative 
correlations,  showing  that  the  London  area  is  not  a  significant  source  of 
PCB's.  The  10-15  km/hr  wind  speed  class  at  station  SAR-S-U1  showed  a  95% 
significant  correlation.  It  should  be  noted  that  while  a  statistically 
significant  correlation  was  found  for  the  NNW  wind  direction  at  SAR-S-U1, 
the  data  are  too  limited  to  indicate  more  than  a  trend,  as  indicated  in 

Figure  2.  Plots  of  SAR-S-I1  data  for  winds  from  the  NW  and  NNW,  and 
SAR-S-U1  for  winds  from  the  NNW  are  given  in  Appendix  5. 

SUD-S-S1 


Winds  from  the  south  gave  a  90%  significant  negative  correlation  at  this 
station.  A  general  positive  trend  with  winds  from  the  northwest  quadrant 
was  found,  but  no  statistically  significant  correlations  were  obtained. 
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SUD-S-R1 

A  95%  significant  correlation  was  found  with  winds  from  the  south-southeast. 
Correlation  for  the  south  and  southeast  are  elevated  but  not  significant 
indicating  the  presence  of  sources  in  a  generally  southerly  direction.  When 
winds  were  from  the  east-northeast  a  90%  level  negative  correlation  was 
found.  Negative  correlations  significant  at  the  90%,  95%  and  90%  levels 
were  found  with  wind  speeds  in  the  0-5,  5-10  and  10-15  km/hr  classes, 
respectively. 

THU-S-R1 


A  significant  pattern  of  correlation  coefficients  can  be  seen  on  the 
histogram  for  this  station.  The  correlation  coefficients  from  the 
south-southwest  direction  through  to  the  north-northeast  are  all  positive 
and  most  are  significant  at  the  90%  level  or  greater.  The  opposite 
directions  show  almost  all  negative  correlation  coefficients.  This  suggests 
that  medium  range  sources  exist  in  a  generally  southerly-to-easterly 
direction.  The  largest  correlation  coefficient  was  0.5284  (sig.  level  90%), 
with  winds  from  the  east-northeast.  Local  sources  d^re  suggested  by  the  95% 
significant  0-5  km/hr  wind  speed  class  correlation. 

THU-S-U1 


The  same  general  trends  found  at  THU-S-R1  were  found  at  this  station  as 
well . 
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T0R-S-S1,  TOR-S-Ul 

Both  stations  exhibited  the  same  pattern  -  generally  positive  correlations 
with  winds  from  the  west  and  north  and  negative  correlation  with  winds  from 
the  east  and  south.  This  suggests  medium  range  sources  to  the  north  and 
west.  T0R-S-U1  showed  a  99%  significance  level  correlation  in  the  0-5  km/hr 
class  and  a  90%  significant  negative  correlation  with  winds  in  the  10-15 
km/hr  wind  speed  class.  This  indicates  trapping  of  PCB's  in  the  city 
envelope  during  periods  of  low  winds.  A  plot  of  T0R-S-S1  data  for  winds 
from  the  north  is  give  in  Appendix  5. 

WIN-S-U1 

A  pattern  of  positive  correlations  with  winds  generally  from  the  west 
suggests  PCB's  are  being  transported  from  the  Detroit  vicinity  into  Canada. 
The  largest  correlation  coefficient  was  for  winds  from  the  west,  and  had  a 
99%  significance.  The  0-5  km/hr  wind  speed  class  showed  a  99%  significant 
correlation.  This  implies  trapping  of  PCB's  in  an  urban  air  mass  over 
Detroit  and  Windsor  during  periods  of  low  winds.  A  plot  of  the  wind  speed 
data  is  jiven  in  Appendix  5. 
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5.0  CONCLUSIONS  AND  RECOMMENDATIONS 

5.1  Conclusions 

Statistically  significant  correlations  of  PCB  concentrations  with  wind 
direction  were  found  for  twenty-two  of  the  twenty-five  PCB  sampling 
stations.  A  total  of  65  significant  wind  direction  correlations  were  found 
out  of  400  wind  direction  correlations  tested.  However,  no  provincial  or 
regional  trends  with  specific  wind  directions  were  found. 

Of  the  correlations  between  wind  speed  class  and  PCB  concentration,  34  were 
founiJ  to  be  statistically  significant  out  of  125  performed. 

Some  of  the  correlation  data  for  wind  speed  and  direction  suggested  the 
presence  of  local  sources  of  PCB's  in  the  vicinity  of  many  of  the  sampling 
stations.  Other  wind  direction  and  speed  correlations  appear  to  indicate 
medium  range  transport  of  PCB's.  No  comment  can  be  made  about  long  range 
PCB  transport  on  the  basis  of  the  correlations  performed  and  discussed  in 
this  report. 

5.2  Recommendations 

As  discussed  in  Section  4.0,  some  of  the  significant  correlations  are  based 
on  limited  data.  As  a  second  survey  was  undertaken  in  June  1980  it  is 
recommended  that  when  this  data  is  available,  it  should  be  combined  with  the 
data  in  this  report  to  form  a  larger  data  base  to  confirm  the  correlations 
already  established. 

It  is  also  recommended  that  the  meteorological  data  be  examined  for  internal 
correlations  by  means  of  multiple  regression  techniques.  In  addition, 
correlations  between  PCB  concentration  and  wind  vector  (made  up  of  direction 
and  speed)  should  be  performed. 

Evaluation  of  large  air  mass  movements  should  be  undertaken  to  assess  the 
potential  for  long  range  transport  of  PCB's. 
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APPENDIX   1 
PCB   DATA 


PCB  DATA 


PCB  STATION:  BUR-S-U1 

DAY   PCB  CONCENTRATION  (ng/m**3) 

1  5.2 

2  10 

3  16 

4  12 
5 

6  1.0 

7  2.6 
8 

9  5.7 

10  2.9 

11  5.4 

12  23 

13  6.4 

14  1.5 

15  5.0 
16 

17  4.7 

18  4.6 

19  5.4 

20  1.2 

21  3.2 

22  2.0 

23  2.1 

24  2.3 

25  2.5 

26  7.1 

27  2.4 

28  3.2 


PCB  DATA 


PCB  STATION:  HAM-S-I1 


DAY     PCB  CONCENTRATION  (ng/m**3) 


1 

16 

2 

15 

3 

15 

4 

7.4 

5 

12 

6 

1.8 

7 

5.8 

8 

5.0 

9 

11 

10 

13 

11 

5.3 

12 

6.3 

13 

9.4 

14 

3.1 

15 

22 

16 

6.3 

17 

3.1 

18 

5.0 

19 

— 

20 

1.1 

21 

1.8 

22 

- 

23 

1.8 

24 

2.6 

25 

1.8 

26 

2.3 

27 

3.9 

28 

2.7 

• 


PCB  DATA 


PCB  STATION:  HAM-S-I2 


DAY    PCB  CONCENTRATION  (ng/m**3) 


1 

- 

2 

13 

3 

17 

4 

18 

5 

- 

6 

5.0 

7 

3.0 

8 

9.0 

9 

4.7 

10 

8.2 

11 

5.0 

12 

11 

13 

4.2 

14 

4.7 

15 

4.6 

16 

4.2 

17 

1.0 

18 

4.6 

19 

- 

20 

3.7 

21 

3.0 

22 

2.7 

23 

9.4 

24 

2.0 

25 

2.1 

26 

3.4 

27 

1.4 

28 

2.6 

PCB  DATA 

PCB  STATION:  HAM-S-U1 

DAY    PCB  CONCENTRATION  (ng/m**3) 

1  9.5 

2  16 

4  12 


6 

7 

9 
10 
11 
12 
13 
14 
15 
16 
17 
18 
19 
20 
21 
22 
23 
24 
25 
26 
27 
28 


1.7 
2.4 
8.8 
2.9 
6.2 
7.6 
9.4 

3.7 

4.0 
1.0 
3.2 
1.6 
2.1 
4.0 
8.8 
5.6 
1.3 
2.6 
2.5 
1.7 
3.1 


PCB  DATA 


PCB  STATIDN:  HAM-S-U2 

DAY   PCB  CONCENTRATION  (ng/m**3) 


1 

10 

2 

13 

3 

1 

4 

9 

5 

- 

6 

3.8 

7 

3.6 

8 

8.2 

9 

18 

10 

5.6 

11 

4.4 

12 

11 

13 

5.0 

14 

7.7 

15 

9.5 

16 

4.2 

17 

5.7 

18 

4.4 

19 

4.1 

20 

1.2 

21 

14 

22 

11 

23 

1.9 

24 

74 

25 

1.2 

26 

3.5 

27 

- 

28 

2.7 

PCB  DATA 


PCB  STATION:  KIN-S-U1 

DAY    PCB  CONCENTRATION  (ng/m**3) 

1  5.6 

2  10 

3  9.5 

4  10 
5 

6  1.1 

7  1.9 
8 

9  4.7 

10  6.0 

11  8.4 

12  6.5 

13  6.3 

14  1.6 

15  8.2 

16  2.1 

17  0.4 
18 

19  5.7 
20 

21  2.4 

22  1.9 

23  4.3 

24  2.4 

25  1.2 

26  3.8 

27  5.2 

28  4.3 


PCB  DATA 

PCB  STATION:  LON-S-U1 

DAY    PCB  CONCENTRATION  (ng/m**3) 
1 


4 


7 


9 

- 

10 

14 

11 

5.9 

12 

- 

13 

3.5 

14 

3.8 

15 

11 

16 

4.7 

17 

6.4 

18 

7.7 

19 

18 

20 

13 

21 

9.3 

22 

- 

23 

11 

24 

6.6 

25 

11 

26 

3.8 

27 

1.5 

28 

0.4 

PCB  DATA 


PCB  STATION:  MIS-S-U1 


DAY   PCB  CONCENTRATION  (ng/m**3) 


1 
2 
3 
4 
5 
6 
7 
8 
9 

10 
11 
12 
13 
14 
15 
16 
17 
18 
19 
20 
21 
22 
23 
24 
25 
26 
27 
28 


11 

14 
15 
8.7 

2.1 
5.8 
4.1 
3.0 
7.2 
4.0 
6.7 
10.2 
2.3 


5.1 
3.5 

1.9 
6.0 
1.6 
1.6 
0.9 
4.3 
3.0 
3.4 
2.1 


PCB  DATA 

PCB  STATION:  MIS-S-U2 

DAY   PCB  CONCENTRATION  (ng/m**3) 


1 

15 

2 

11 

3 

10 

4 

5.3 

5 

- 

6 

1.1 

7 

3.4 

8 

7.0 

9 

3.5 

10 

10 

11 

6.9 

12 

5.5 

13 

10 

14 

0.8 

15 

- 

16 

7.3 

17 

1.4 

18 

3.2 

19 

4.5 

20 

- 

21 

1.8 

22 

0.8 

23 

1.7 

24 

2.8 

25 

1.4 

26 

2.0 

27 

4.4 

28 

1.3 

PCB  DATA 


PCB  STATION:  MIS-S-U3 


DAY   PCB  CONCENTRATION  (ng/m**3) 


11 

11 

18 

11 

15 
5.6 
2.7 
9.0 
2.7 

21 
5.8 
8.7 
9.9 
3.2 


4 
5 
6 
7 
8 
9 

10 
11 
12 
13 
14 
15 
16 
17 
18 
19 
20 
21 
22 
23 
24 
25 
26 
27 
28 


2.3 
3.0 

2.2 
2.5 
0.6 
1.6 
0.8 
2.7 
2.1 
2.2 
1.6 


PCB  DATA 

PCB  STATION:  MIS-S-U4 

DAY    PCI3  CONCENTRATION  (ng/m**3) 


1 

12 

2 

7.2 

3 

9.4 

4 

3.1 

5 

- 

6 

0.9 

7 

5.0 

8 

5.3 

9 

2.6 

10 

11 

11 

7.7 

12 

15 

13 

8.9 

14 

3.5 

15 

5.5 

16 

5.0 

17 

4.0 

18 

4.5 

19 

2.3 

20 

9.0 

21 

4.7 

22 

1.3 

23 

1.6 

24 

3.3 

25 

3.4 

26 

3.3 

27 

3.7 

28 

1.4 

PCB  DATA 

PCB  STATION:  MOR-S-U1 

DAY     PCB  CONCENTRATION  (ng/m**3) 


1 

16 

2 

6.8 

3 

2.6 

4 

3.9 

5 

12 

6 

6.0 

7 

1.0 

8 

1.8 

9 

4.8 

10 

5.2 

11 

6.4 

12 

8.6 

13 

7.3 

14 

0.6 

15 

- 

16 

6.4 

17 

8.6 

18 

6.8 

19 

1.8 

20 

2.2 

21 

1.8 

22 

3.2 

23 

- 

24 

3.0 

25 

- 

26 

2.3 

27 

1.8 

28 

3.9 

PCB  DATA 

PCB  STATION:  NAN-S-R1 

DAY    PCB  CONCENTRATION  (ng/m**3; 


1 

7.5 

2 

9.3 

3 

11 

4 

- 

5 

14 

6 

9.5 

7 

4.3 

8 

12 

9 

9.6 

10 

7.8 

11 

3.0 

12 

5.0 

13 

7.7 

14 

1.0 

15 

3.4 

16 

4.8 

17 

2.0 

18 

4.9 

19 

1.6 

20 

1.6 

21 

2.2 

22 

1.1 

23 

1.5 

24 

2.4 

25 

3.3 

26 

3.2 

27 

3.6 

28 

2.5 

PCB  DATA 

PCB  STATION:  NAN-S-R2 

DAY    PCB  CONCENTRATION  (ng/m**3) 

1  14 

2  9.2 

3  12 

4  9.4 

5  18 

6  1.8 

7  1.3 

8  4.9 

9  2.3 

10  7.1 

11  4.4 

12  7.9 

13  4.7 

14  1.1 

15  2.5 

16  7.8 

17  2.6 

18  5.4 
19 

20  2.5 

21  2.1 

22  3.9 

23  2.1 

24  3.8 

25  2.3 

26  1.8 

27  0.9 

28  2.4 


PCB  DATA 

PCB  STATION:  OSH-S-U1 

DAY   PCB  CONCENTRATION  (ng/m**3) 


1 

6.6 

2 

14 

3 

12 

4 

10 

5 

- 

6 

1.2 

7 

4.1 

8 

7.1 

9 

14 

10 

12 

11 

7.3 

12 

10 

13 

12 

14 

1.2 

15 

— 

16 

3.5 

17 

1.0 

18 

2.2 

19 

4.8 

20 

1.2 

21 

2.2 

22 

1.7 

23 

3.5 

24 

1.0 

25 

6.3 

26 

3.5 

27 

2.6 

28 

2.0 

PCB  DATA 

PCB  STATION:  SAR-S-I1 

DAY     PCB  CONCENTRATION  (ng/m**3) 

1  20 
2 

3  21 

4  7.8 

5  17 

6  10 

7  3.7 

8  10 

9  6.6 

10  6.5 

11  3.8 

12  11 

13  9.8 

14  0.6 

15  3.5 

16  8.3 

17  4.4 

18  2.0 

19  2.6 
20 

21  2.2 

22  4.9 

23  3.0 
24 

25  6.5 

26  8.3 

27  1.1 
28 


- 


PCB  DATA 

PCB  STATION:  SAR-S-U1 

DAY    PCB  CONCENTRATION  (ng/m**3) 


1 

- 

2 

17 

3 

47 

4 

- 

5 

- 

6 

7.4 

7 

5.0 

8 

10 

9 

6.2 

10 

8.6 

11 

5.7 

12 

26 

13 

13 

14 

5.2 

15 

- 

16 

9.0 

17 

4.1 

18 

12 

19 

7.8 

20 

5.0 

21 

1.8 

22 

4.6 

23 

8.0 

24 

9.0 

25 

3.8 

26 

4.4 

27 

9.3 

28 

6.2 

PCB  DATA 

■ 

PCB  STATION:  STC-S-R1 

DAY     PCB  CONCENTRATION  (ng/m**3) 


1 

19 

2 

8.3 

3 

23 

10 


6 

2.2 

7 

2.1 

8 

- 

9 

5.7 

10 

12 

11 

4.7 

12 

6.4 

13 

6.9 

14 

5.9 

15 

8.0 

16 

4.4 

17 

1.9 

18 

• 

19 

2.0 

20 

3.2 

21 

11 

22 

1.6 

23 

3.2 

24 

1.8 

25 

2.6 

26 

2.9 

27 

2.4 

28 

4.3 

PCB  DATA 


PCB  STATION:  SUD-S-S1 


DAY    PCB  CONCENTRATION  (ng/m**3) 


1 

4.3 

2 

12 

3 

9.2 

4 

5.0 

5 

- 

6 

4.2 

7 

1.1 

8 

6.7 

9 

- 

10 

8.1 

11 

6.3 

12 

5.2 

13 

9.4 

14 

3.7 

15 

- 

16 

- 

17 

2.2 

18 

5.0 

19 

2.3 

20 

2.4 

21 

2.6 

22 

1.0 

23 

1-4 

24 

2.1 

25 

1.5 

26 

3.2 

27 

2.3 

28 

3.9 

PCB  DATA 


PCB  STATION:  SUD-S-R1 


DAY    PCB  CONCENTRATION  (ng/m**3) 


1 

15 

2 

12 

3 

9.6 

4 

3.9 

5 

7.7 

6 

4.2 

7 

2.9 

8 

6.2 

9 

- 

10 

5.6 

11 

11 

12 

5.9 

13 

10 

14 

3.0 

15 

3.5 

16 

1.8 

17 

- 

18 

2.6 

19 

1.9 

20 

2.3 

21 

6.4 

22 

1.7 

23 

8.7 

24 

16 

25 

10 

26 

3.0 

27 

3.0 

28 

17 

PCB  DATA 

PCB  STATION:  THU-S-R1 

DAY    PCB  CONCENTRATION  (ng/m**3) 


1 

8.8 

2 

9.3 

3 

7.1 

4 

2.0 

5 

- 

6 

3.5 

7 

0.7 

8 

- 

9 

3.3 

10 

4.3 

11 

4.8 

12 

10 

13 

4.7 

14 

2.5 

15 

2.5 

16 

2.4 

17 

5.8 

18 

4 

19 

2.5 

20 

2.2 

21 

2.5 

22 

3.6 

23 

1.3 

24 

0.8 

25 

1.0 

26 

1.9 

27 

2.4 

28 

3.2 

PCB  DATA 


PCB  STATION:  THU-S-U1 


DAY 


PCB  CONCENTRATION  (ng/m**3) 


3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

23 

24 

25 

26 

27 

28 


13 

7.7 

14 

12 
9.9 
5.0 
1.6 

13 

11 
3.9 
5.7 
6.8 
7.9 
3.7 

1.7 
11 
6.8 
1.7 
1.7 
1.5 
1.3 
2.2 
2.7 
2.6 
1.3 
1.5 
4.3 


• 


PCB  DATA 

PCB  STATION:  TOR-S-S1 

DAY    PCB  CONCENTRATION  (ng/m**3) 

1  11 

2  10 
3 

4  15 
5 

6  2.2 

7  3.8 

8  5.8 

9  4.6 

10  8.2 

11  4.0 

12  6.3 

13  7.3 

14  2.9 
15 

16 

17  3.1 

18  3.5 
19 

20  1.9 

21  5.1 

22  4.6 

23  2.6 

24  3.5 

25  5.0 

26  4.5 

27  4.4 

28  1.6 


PCB  DATA 


PCB  STATION:  TOR-S-U1 


DAY 


PCB  CONCENTRATION  (ng/m**3) 


1 
2 
3 
4 
5 
6 
7 
8 
9 

10 
11 
12 
13 
14 
15 
16 
17 
18 
19 
20 
21 
22 
23 
24 
25 
26 
27 
28 


14 


11 
25 

5.9 


16 
7.9 
3.0 
5.9 

11 
5.2 
7.8 


5.6 

2.3 
14 
3.4 
4.5 
1.9 
6.3 
6.0 
3.9 
5.2 


PCB  DATA 

PCB  STATION:  WIN-S-U1 

DAY    PCB  CONCENTRATION  (ng/m**3 


1 

6.6 

2 

11 

3 

13 

4 

10 

5 

16 

6 

4.8 

7 

2.2 

8 

5.9 

9 

9.5 

10 

6.1 

11 

8.5 

12 

5.2 

13 

6.0 

14 

5.6 

15 

6.8 

16 

2.9 

17 

3.2 

18 

4.5 

19 

2.8 

20 

2.7 

21 

4.3 

22 

2.0 

23 

- 

24 

2.2 

25 

1.5 

26 

1.5 

27 

3.3 

28 

3.4 

APPENDIX  2 
COMPUTER  PROGRAM  LISTING 


"^.  1 


,-  no  r  ir  |_  at-  i-        p  -£: 


K-r       -T->  .  -      K'—  -      r.£T_5 


"OR    THE    2£ 


-::-    WIND       SF'EEi       DATA 


■ijrrQ     i|CC      T  C;     '"^  f\S£rf     p 


OR  CHANCE  TL'.Z 


T  M  C ' 


-at!  w  -  —  v;C'C- 


■J  rr 


/.  PnTKJT:-;  OF  CCMPA 


INS   DIRECTION? 


•"'.   SPDCLASS    -  CLASSIFIES  WINS  SPEED  INTO  5 

H.  CORRELATION  -  CALCULATE'?  CORRELATION 
COEFFICIENT  AND  SIGNIFICANCE 

I.  3UTPU7  -  OUTPUTS  CAi  CL'LATED  CORRELATION 
COEFFICIENTS  AND  ASSOCIATED  STATISTICS 

J,  GRAPH  -  OUTPUTS  HISTOGRAM  PLOT  OF 
CORRELATION  COEFFICIENTS  VS  WIND 
DIRECTION  CLASS 

K,  CLOSE 'FILE  -  STCRES  CURRENT  PCE  OR 
STATION  DATA  ON  A  DISK  FILE 


1-.       :     i"\  '_ 


—        .-  C 


*CHES    PCS    AND/Cr    STATIC* 


A    DIOr.    FILE 


•/,'P  «     l.'MT  £  '  *  *?) 


-,  -,      O  /'  1 


o  ae    <- 


"i.    ::n" 


V i  t*P-"  i    '  i  *  *• 


»       ri  /an 


PK'P  '• 


MAS  ■) 


M    PC2NAM*  (25 )  *  tfrTNAM*  (13: 
M    TL  I  NEW  INDC  (16)?  St-.'  I NDD  C  1 
M    PCS  (28)*  STATION* 29,  24-  3 
"-:    VC23}*  DIRCLAS(28,  16),  3P 


-  3 1 GND IP* (  16: r  SZ SKSFD 5(5; 

_  ■  r  TLIK'EWIND"  '  2:  »  SWINDS  (  51 

) ,  AVTE;*P(2S3  r  DEVTEMP<2e))  *  X  (29) 

DCLAS(2S,  5 '7  1,'DCOEFFC  16) 


nA~a  c  h 


L 


E-._   GRA: 


-•.-  E 


rH*  (  32?  1  )  *  GRA  - 'H*  ( 32*  4  )  *  GRAPH*  ( 32r  5 )  * 
K$C32*  10)*  £RAPH*{32:  14) '*  GRAPH*  (32*  i 
.jr'Hf  ( 32*  24 ) ,  2RAPH*  (32.-  25)  *  GRAPH*  (32*  2 

PEAT'  GRAr>H*(32*  50),  GRAPH*  •'-2*  51)*  GRAPH*  (32*  3 


p— &n 


GRAPH*  (32*  60),  GRAPi 


P?iPH4  '  ■""' 


(j«t  ( 


:-?.k"H*(32»  74)  *  "2 


GRAPH* ( 32:  S ) *  GRAPH* ( 32*  9 ) 

5; *  GRAPH* (32*  16;* GRAPH* (22.  20; 

£ ) ,   GRAPHS (32r  2?) *  GRAPH* ( 32.-  30) 

' ),  GRAPH* (32* 55) *  GRAPH* (32* 5?) 

3  J ,  GRAPH*  C  32*  £■?  >  *  GRAPH*  (  32*  70 ) 


DATA    7HU-  :--.!i  7Hl> S-R 

DATA    .--'A.|wi--S-  I2f  :-'AM-2-'Jl,  HA^-f 


DATA    K0E-77C23*  ASS-    TE-AY,  KYD- 


2034*  AEG-    LOHDt  ^L'E-29026 


DATA 

DA  FA       GPGGLAGGr  OGRRELAriC!".'-  31 


"PUT*  GRAPH*  CIGGE'FILE,  RETRIEVE 


T  =  "i      Ti"      "  — 


.-A 


REA2    "!ETMA^r  (  I ) 

T 


.::  «:,  ff stats,  ga 


P A] M  —<;"'*#!        *    '-   .    f~  : '•.''■'■• 


*  -I 


"AD  ( -  1  •  0 ) 


-  — v  I 

END 


-  q  r, ;  ( - 1  .  o  i 


ECHO,  A\fAi  YZEt  fiddle: 


SUGRSUTINZ  DATA ' ENTER 

w.  Jj.  <;.  v.  g,  ;<  K  v.  l:-  ^ . .-    ;-  ■-;.  :_-  :-  -x.  *.  .*.  :  •>  ■';-  <-  ■>  -i-  ';-  -i  #  <-  *•  X-  ■«■  >  >-  >#■■£■■*■■£■•»  «■  £  >-  >  -"-  ■»*  #■  #  *  i-  :-'r  V 


ER!  INPUT    "ENTEK    ~C3     (1)     DP     E-TATION    C2.'     DATA*     OR    QUI"      ( 3  5  %  A 1 

T~    ^  ■'=:":    R— T' 'RN! 

DM    A!     30TG    PCB'EN'TRr  STAT 'ET.TR 
GOTO    DATA 'ENTER 

IMF'U  ;     " 'r'TA  FI  L'i\    ':;TART     >  .ME*        >  &  i  AR  . 
r:-    1  =  1    TC    tfAX'/SATKER-l 

:kput  pcs a  3 

'."~VT  T 

CALL  CLOSE 'FILE 
30TC  DATA 'ENTER 

'ENTR:  INPUT  "NUMBER  CF  KEU    ".E7    DATA  FILE  M  -  IS?'?  "tCTATNUM 
INfPUT  "DAY  ~C  START  ENTER 1X3  TEMP  DATA:  "fSTARTMY 
:F  ST ARTDAY>M.*X WEATHER  OCTO  IN^I.VD 

IF  STARTDAYti  FLA'32=3:  A2=3TATNL'M:  I=STARTDAY:  CALL  RETRIEVE 
70R    I=STARTDAY  TC  MAXWEATh'ER 
FOR  .J=l  TO    24 

PRINT  USING  "DAY  ##   -  ", I J 

SOT  KIT     II-TV~.      wphfTSTO     TMOCCj" — '  re  —     CftC1     P«"»J  «C      *•£  ji  •       "  .  ,  1—1; 

INPUT    "    ",  STATION* I*  J»  1  J 

STATION  (I,  J,  1  > -STATION*  It  J,  D/IO. 
NEXT    J 

P R I T ,'  T    CK R  t  (  0  0  7 ) 
NEXT    I 

CALL  CLOSE7 FILE 

IF  =:TARTDAY#:  :::?U7    "'DAY  TC  START  ENTERING  SPEEE'  DATA:  %  S 
IF    STAR  rDAY  >M  AX!-.1  EAT  HER  "ITS  DIR3N 
-p  5~Y.~~~AY#!  FLAG2=2" A*I=.STATfsL'!*.,rCAi.L  RE~R' IE  'E 
FOR  I=STARTDAY  TO    ?**  AX  WEATHER 
FOR  J=l  TO  2 A 

PRINT  USING  "DAY  ##  -  :r,  I, 

PRINT  USING  "EXTEP    WIND  IRZED  FOR  HOUR  ##:  ",  ■_'-:? 
INPUT  "  %  STATION*  I,  J*  35 
NEXT  J 

NEXT  I 


:RY 


Fi  I  s    Z  I '. ..  1 :  PC?.  BAS  l  2<~.  ■"  r  '  1    pfi  n  te 


' '-  •  > 


riSS     4 


■—( 

-.;_" 'U  —  2 

.?■ 

_L     ?lc:i  '.- 1 

F 

S7APTI  Wv  -• 

— 

,-  —  ,-.— -.  y  -.    M 

r.xwEr 


"DAY 


-J-l 


skipdirn: 


ME  XT    ..' 

r,tA  I 


JH  !  h  '  S.IM 


mtf; 


i  h,  £.  ^_ .-».  v-  >  -.-  ■£ 


••_  .--  -  .  >'   .-.  -X.  M,  .X.  -'-.  &.  -V. 


^C>..;     -;-    fc.^C.M.^.    ;-.X.    ###Jfc#fc#    JHR?*r#-8-**-fc4&4S#£r*i### 


IF  -~2  RETURN 
INPUT  "DO  YOU  WAV 

if   Ai^:iri:="v 


b*t->1 TV   I   -     '■ 


DATA, 


I  :.iKf'ATA: 


1  •- 


wnr  rev 


FOB ' MOD 


'AT 'MOD 


rifiV'-  "iSTfl  nn  v  ~ ' .'  WANT  ri~!  MODIFY?  '%  DAY 
ON  A 1  GOTO  FOS '  M 3B»  STAT  '  K  3D 

INPUT  "ENTER  NEW  POD  CONCENTRATION:  ">  PCBCPAY3 
GOTO  DATA ' SAME 

INPUT  "MODIFY  "If-'D  (i>»  DIRECTION  COS  OF  SPEED  (3)?  "r  A3 

INPUT  ,r7:u:E  T'Z    FE  MODIFIED"  '%  TIMEDAY:  TIKEDAY-TIMEDAY+i 

FPINT  USING  "CURRENT  VALUE  IS:  £*'#*'.  #r,r  STATION*  CIAYr  TINZDAYf  A33 

INPUT  "ENTEF  THE  NEW  DATA: " , STATION ( DAVr  FIMEDAY, A3 ) 
3CT0  DATA '"-AVE 

INPUT  "SAVE  THE    DATA  ON  DISK?  %  A 2* 

IF  A~S=:"Y"  OF  A2S=,rYES"  CALL  CLOSE 'FILE 

3070  CORRECT 


■fc-  :r  -:>  --  >-  ----/-  -  -     >  ■>■    >  <-  -*•  *L  #  ---  '-.   '-'  *-  '•-  ' 


,:  :..-.  :      -.  3     ;.  '.--•    >    «fc£    :...     -.-  >     .-  ..  ,-  >.   ;-:-•>-  X-  <  ■>  >---:-*■  *•  -:-  i  £  -:-  '-  >-  <  '•  * 


•  v  %  #  X* 


y,        ■   .  <   :.  >.  «.£.#*  *  *-  = 


i  !v 


•~P 


iiHi -  _ 


C23     DATA,     OR    SUIT    (S3     ".-Ai 


• 


».■   i  j.     r  :r 


K-  T  —  J> 

.  i  •  -    . 

'•     :"  1 


Hi   '  i_ . 


'r'CZ    STATION    NUNE'EP?    (1    -    26): 


l-' '  "   IV 


:all  retrieve 

:HOPC?:    PRINT    CHR$(0i23 

;---    i  =  :    re   ?: print   "    "smexT   I 
~"::::~  us  INC-   ";  fr    -^-t-v 


T        " 


C'rOVAM 


PCBKUH) 


rvin 


:  ps i is  i 


or  :=i    re  wAxwEArKEK-i 

PRINT   USIN3    "  #E#Mr  I J 

-jr    PC2M)£9*    PSI N!7    '-'rZNG    F3l"rTAB(l 


-  _.  i  -     .-.ma 


TC»T         "-UOA   i"   •■•■  1    "t 


CALL  retrieve 

hT:FCR    1=1    TO    MAX WEATHER 
PRINT    CHR$<012) 

for  n=:  ra  3: print  ,!   "::next  :: 

PRINT  FA3  (  31  >?  "MET  STATION  DATA" 


^:r  :i  =  i  to  4: print  !i  ":next  is 

r  .  .  i  r-i  .   _•  r-  j.  r\L  up.  :  ?l ..-  ;f  f  x 

"RlKfT  T&Si  13)  ?  "  CDEG  C;  :'  ?  TAB  (  23)  ;  "  'KN/HR  )  a  •  TA3C37)  ?  & 
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PRINT    7AS< 19) ?  'TORRE 

3RAPH$  f  !»•  I  >~WI.N 
CRAPH?  (2.-  I  )=PC3 
NEXT    : 
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-  |  r  -  t  ,  - 1 2 
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IF    G  =  r:-iGT^0    THEN    GO    TO    SKIP2: 

Trr     i^pUj-inj-rwj     OFKKGT^GFHKG"'";- 1 
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INPUT    "STATION   FILE!    NUMBEF 
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CONTINUE-  .  .  "-  DUMMY 
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APPENDIX  5 
PLOTS  OF  SIGNIFICANT  CORRELATIONS 
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